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The initial value problem of scalar-tensor theories of gravity (STT) is analyzed in the physical 
(Jordan) frame using a 3+1 decomposition of spacetime. A first order strongly hyperbolic system 
is obtained for which the well posedness of the Cauchy problem can be established. We provide 
two simple applications of the 3+1 system of equations: one for static and spherically symmetric 
spacetimes which allows the construction of unstable initial data (compact objects) for which a 
further black hole formation and scalar gravitational wave emission can be analyzed, and another 
\ application is for homogeneous and isotropic spacetimes that permits to study the dynamics of the 

Universe in the framework of STT. 
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I. INTRODUCTION 



A physical meaningful theory is supposed to have a well posed Cauchy problem. Otherwise its predicting power 
during the Cauchy development can be lost. 
■ In general it is not a trivial task to determine if a field theory possesses a well posed initial value problem. The 
systematic approach is to write the field equations as a manifestly hyperbolic (e.g. symmetric or strongly hyperbolic) 
^sO . system of partial differential equations (PDE). The existence and uniqueness theorems of such a system guarantee 
(under quite general hypothesis) that the solutions exist are unique and depend continuously on the initial data 1 . 
In the case of gauge theories the analysis of the Cauchy problem becomes rather subtle or even complicated due to 
CN ■ the gauge freedom of the theory. Indeed, the hyperbolicity of the PDE associated with gauge theories usually depends 
' on the gauge choice. One must ensure that the chosen gauge is well behaved in the sense that it is preserved during 
the evolution if it holds in the initial spacelikc hypcrsurfacc. 

As regards general relativity (GR), it took several years and considerable effort to establish the well posedness of 
the Cauchy problem in the so-called harmonic gauge (e.g. see Refs. |lj). In fact this gauge allows to write the Einstein 
field equations as a set of quasilinear diagonal second order hyperbolic system for the metric components (the reduced 
field equations). The Bianchi identities ensure that the harmonic gauge is preserved during the Cauchy development 
if it holds initially. 

During the last three decades a considerable work has been made for solving the Einstein's equations numerically. 
Unfortunately, several experiments have revealed that the harmonic gauge is not the best coordinate choice to evolve 
initial data that culminates in the formation of black holes. 

On the other hand, during the past few years, several novel formulations of Einstein's equations have been proposed 
for treating the Cauchy problem numerically. Unlike the traditional ones (notably the ADM formulation) these new 
formulations are strongly hyperbolic (see Ref. @, Q for a review). This analytic improvement together with new 
gauge conditions have dramatically enhanced the numerical stability during the evolution of several initial data. 

From the observational point of view there is little doubt about the validity of GR. So far almost all the predictions of 
the theory have been corroborated by experiments. It is expected that during the next decade or less the gravitational 
waves predicted by GR be directly detected by several interferometers. This information can be crucial for further 
testing the theory in the future. 

Now, despite that GR seems to be the simplest diffeomorphic gravitational theory and so far the best one as 
observations are concerned, there are some astrophysical and cosmological data whose interpretation has renewed the 
interest for alternative theories of gravity. In the cosmological setting, the discovery of the accelerated expansion of 
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the Universe together with the recent probes of the cosmic microwave background, points to the existence of some 
unknown type of matter (termed dark energy) which exerts negative pressure. The simplest implementation of this 
dark energy is to introduce a cosmological constant in the field equations. However, for several reasons (for instance, 
the so called coincidence problem) this is not considered to be the best solution to this enigma, and therefore cosmic 
dynamical scalar fields (termed generically as quintessence) have been proposed as alternatives for explaining such an 
acceleration with the property that their energy-density evolves with the scale factor. 

Very recently it has been advocated that modified theories of gravity can produce a cosmic acceleration without the 
introduction of a cosmological constant or of any other form of dark energy [3| • These theories seem however to face 
serious problems at the solar system level [|| and show different kinds of instabilities [f| 0] ■ 

At smaller scales (galactic scales) the dynamics of stars and gas around spiral galaxies shows that their tangential 
velocity around the bulk of the galaxy is independent of the distance to the center. These observations contradict 
the Newtonian dynamics that was supposed to be valid at such regime. The explanation of this data (termed as 
flat rotation curves of galaxies, hereafter FRC) remains a mystery. The standard approach to solve this problem (as 
in cosmology) is to assume that another form of matter (dark matter) generates a galactic halo that interacts only 
gravitationally with the visible matter. This dark matter would produce an energy-density which decays with distance 
as l/r a (with a w 2) such that the mass of the halo grows more less linearly with r and therefore implies (within a 
Newtonian framework) a tangential velocity which is almost constant. 

There exists also the view that the introduction of such dark matter does not really solve the problem in that the 
nature of such matter is still unknown. In this direction, the most popular and successful alternative explanation to the 
FRC without the need of any dark matter is the so-called modified Newtonian dynamics or MOND L S] , which proposes 
that Newtonian dynamics is only valid for "large" accelerations while for small accelerations the Newtonian gravity 
is compensated by a dynamics which is proportional to quadratic accelerations. Therefore, in the low acceleration 
regime, the modifications give rise to tangential velocities that are independent of the distance. 

For many years the lack of a relativistic generalization for MOND was one of the main objections by its detractors. 
Very recently, Bekenstein suggested a tensor-vector-scalar theory of gravity (TeVeS) [9j which is generally covariant 
and which reduces to the MOND model in some limit. Needless to say, this theory should be considered as a viable 
alternative theory as long as it is compatible with all the observations (cosmology, solar system, binary pulsar, etc.) 
not only those related with the FRC. On the other hand, recent analyses indicate that TeVeS theory exhibits certain 
kind of instabilities too 0] (see also Ref. [l(| for several criticisms to TeVeS). 

As mentioned in the beginning of this introduction, a serious candidate to replace GR should also have a well posed 
Cauchy problem. Therefore, in addition to the confrontation with observations, this is a task that should be tackled 
soon if some of these alternative theories are to be considered as serious theories. 

Scalar-tensor theories of gravity (STT) arc another alternative theories. The most prominent example of STT 
is perhaps the Brans-Dicke (BD) theory STT introduces a new fundamental scalar field which appears to be 

coupled non- minimally (NMC) to gravity (in the so-called Jordan frame). Due to this feature STT predict several 
phenomena not present in GR. STT are perhaps the best analyzed alternative theories of gravity. When the effective 
BD parameter is sufficiently large, the STT are compatible with the solar system tests. Moreover, STT have been 
proposed for explaining the acceleration of the Universe without any cosmological constant fl2l. Ha . fl4j| . 

Finally, from the mathematical point of view, STT posses a well posed Cauchy problem [la . 1 1 €31 ] . 

One of the main motivations of this contribution is to further discuss the Cauchy problem of STT in the Jordan 
frame [l6| . In fact, it is customary to find in the literature remarks about the unsuitability of the Jordan frame to 
analyze the Cauchy problem [l7|. This widespread prejudice is often invoked to transform the action (TT|) (see the 
next section) to the so-called conformal or Einstein frame where the NMC function F(<fi) is absorbed in the curvature 
by defining a non-physical conformal metric (the metric which is not associated to the non- null geodesies). The 
mathematical advantage of the conformal frame is that the resulting field equations, notably the equations for the 
non-physical metric, contain an effective energy-momentum tensor which does not include second order derivatives of 
the scalar field. This feature allows one to analyze the Cauchy problem of STT in the Einstein frame more or less 
in the same way as in GR, since basically the first order derivatives of the scalar-field which appears as "sources" 
do not contribute to the principal part of the system of equations. Of course the price to be payed when working 
within the conformal frame is that one is always dealing with non-physical fields (in the above sense) and that the 
transformation back to the physical frame has to be well defined. This is not always the case, in particular when F((f>) 
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is not positive definite one can encounter situations where the conformal transformation is singular [1 81 ] 2 . 

In this contribution we shall first remind the field equations of STT treated in 3+1 form Then we sketch the 
hyperbolicity analysis that leads to a well posed first order system of PDE. 

In section 4, we consider the very simple case of initial data that represents a static and spherically symmetric 
spacetime. When such data correspond to an unstable configuration (e.g., a neutron star with maximal mass), the 
data can be used as a departure point for studying the gravitational collapse of compact objects within STT and the 
emission of scalar gravitational waves. A novel system of equations is presented using area-r coordinates (spherical 
coordinates where the area of 2-spheres is given by 4irr 2 ) and the system reduces to the one usually employed in GR 
in the absence of a scalar field. 

We provide in section 5, the fundamental equations required to construct an isotropic and homogeneous cosmological 
model for the Universe in the framework of STT. 



II. SCALAR-TENSOR THEORIES OF GRAVITY 



The general action for STT with a single scalar field is given by 



S[g ab ,<p,tp] 



16ttG 



R 



i(v^) 2 + u(</>) 



-gd 4 x + S ma ,tt[g a b,ip] 



(1) 



where -0 represents collectively the matter fields (fields other than cj>; units where c = 1 are employed) 3 . 

The representation of STT given by Eq. JT]) is called the Jordan frame. The covariant field equations obtained 
from variations of the action with respect to g a b and <fi are: 4 
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(5) 
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Using Eq. (J5]), the Ricci scalar can be expressed in terms of the effective energy- momentum tensor (EMT) Eq. Q 
and then Eq. takes the following form, 



/U'-2/'U-i/'(l + 3/")(V^) 2 + i/'T matt 
U< P - 7 T77I\ ' i 8 ) 

f{ 1 + 3 4r) 



2 In a very recent analysis [XSj] , the authors argue that the connection between both frames also breaks in the limit of GR due to a lack 
of differentiability of the scalar-field redefinition in this limit. In fact the authors show that the cosmological attractor property of GR 
with respect to STT [2q | depends on the choice of the frame. 

3 Although the function F(<f>) is, in principle, completely arbitrary, it turns from the analysis of the Cauchy problem, that F(cj>) > is 
a necessary condition to avoid potential singularities or blow-ups in the equations and in several variables (this condition corresponds 
later to /(</>) > in the notation adopted in this section). 

4 Latin indices from the first letters of the alphabet a, b, c, ... are four-dimensional and run — 3. Latin indices starting from i k, ...) 
are three-dimensional and run 1 — 3. 
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where T mat t stands for the trace of T^ att and the subscript "matt" refers to the matter fields (fields other that 
Now, the Bianchi identities imply 

V c T ca = . 

However, the use of the field equations leads to the conservation of the EMT of matter alone 

V T ca — 



(9) 



(10) 



which implies the fulfillment of the (weak) equivalence principle. 

In which follows we shall use the 3+1 formalism of spacetime [U in order to recast the field equations as a 
Cauchy initial data problem. 

The definition of the 3+1 variables and the evolution equation obtained from these are respectively: 
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where Cn is the Lie derivative along n which is the unit time-like vector field normal to the spacelike hypersurfaces 
St used to foliate the spacetime (as long as it is globally hyperbolic); 7 is the 3-metric induced on E t and D is the 
covariant derivative compatible with the 3-metric. In components n a = (1/ a, — f3 l / a) where a and (3 1 are the lapse 
function and the shift vector, respectively, associated with the spacetime coordinates (t^x 1 ). 

Moreover, the 3+1 formalism allows one to write the field equations as initial data constraints and first order in 
time evolution equations (see Ref. [l6j]) 5 . The field equations ([2]) then read as the following constraints (H~5l) and (fTB]) 
and the evolution equation (fTT|) : 
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where £ matt := n a n b T^\ S^ tt := J a c 7 b d T^ att , S matt := S a a matt , J£ att := -n dl a A t and 7 a c = S a c + n a n c is 
the projector onto S t . 



5 A slightly different notation with respect to [ill has been employed here in order to match with the notation used in many references 
on numerical relativity. In order to return to the notation of Ref. [JJ|, one needs to perform the transformation a — > N, f) 1 — + — N l , 
7ij — > h%j for the lapse, shift and the 3-metric respectively. 
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From Eqs. (fT5|) and (fTT|) , a useful evolution equation for the trace of the extrinsic curvature can be obtained: 
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where 3 A is the Laplacian compatible with the 3-metric 7y, and Q 2 = QiQ 1 . 
The evolution equation for II is obtained from Eq. © and reads as follows 16], 
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matt — *^matt -^rnatt 



III. HYPERBOLIC SYSTEM OF EQUATIONS 

The "gravitational" sector of the equations (fTS)) — (fT51) (i.e., the parts that remain in the absence of the scalar field) 
corresponds to the usual ADM equations a la York [2JJ (hereafter AD MY system). These equations do not form 
a strongly hyperbolic system, but only a weakly hyperbolic one, for which the Cauchy problem is not well posed . 
Therefore, the inclusion of the non-minimally coupled scalar field will not improve the situation. 

In pure GR various strongly hyperbolic systems can be constructed from the AD MY equations [3J]. The approach 
consists in recasting all the equations in full first order form (in space and time), and then to add appropriately 
multiples of the constraints to the evolution equations. In this way one can obtain a system in the form 

d t u + M l V l u = S{u) , (20) 

where u represents collectively the fundamental variables (like the /ly's, -fiQj's, etc.), M 1 is called the characteristic 
matrix of the system (which depends on u but is independent of its derivatives) and S(u) are source terms which 
include only the fundamental variables (not their derivatives). 

The addition of the constraints can make the characteristic matrix to be diagonalizable (i.e., the matrix has a real 
set of eigenvalues and a complete set of eigenvectors). 

Although there are linear-algebraic methods (see for instance [23|) to prove the strongly hyperbolic character of the 
characteristic matrix, one can also follow another economical but equivalent approach. This consists in constructing 
the eigenfields w := M _1 m from the eigenvector matrix M with the property L = IR^IVPIR. is the diagonal eigenvalue 
matrix associated with a particular direction l x' of propagation. The point is that one can achieve such a construction 
by a judicious guessing method (see the Appendix for a simple example). This is roughly speaking the approach 
followed by Bona-Masso [2~jj . 

Of course, another key element in the construction of a strongly hyperbolic system is the gauge choice. In our case 
we consider the shift as a given function of the coordinates and a modified Bona-Masso time slicing which reads 

d a = -h M {a)a 2 [K - -^-^II) , (21) 

\ JBM / / 
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where do — dt — /3 l di, /bm(«) > (i.e., /bm is a positive but otherwise arbitrary function of the lapse) and is a 
parameter 6 . 

In pure GR the usual BM condition [i.e., Eq. (|2"Tj) with 9 = 0] generalizes the harmonic slicing ]bm = 1- The 
BM condition therefore includes several choices of /bm(«) that can be more useful than others depending on the 
evolution to be analyzed (see |3| for a more detailed discussion). However, the usual BM condition seems not to work 
for STT in the Jordan frame. In fact it is possible to show [25[ that taking = leads to a weakly hyperbolic system. 
Therefore, the modification of the usual BM slicing condition is a crucial step towards a strongly hyperbolic system 
of STT when formulated in the Jordan frame. 

In this way, the evolution system in terms of the eigenfields and their corresponding eigenvalues reads [25j 

d Q w l + d x (AW) ~ (22) 

(here no sum over the repeated index) where the symbol ~ indicates equal up to principal part; w 1 represents 34 
eigenfields constructed upon the first order variables defined in Sec. II, and also from first order variables defined in 
terms of the spatial derivatives of the 3-metric and the spatial derivatives of the lapse. In pure GR there are only 
30 eigenfields (see Ref. [|| for a review) since Qi = = II. Here A' is the eigenvalue which is associated with the 
speed of propagation of the corresponding eigenfield vf . In fact there are 20 eigenfields which evolve along the time 
lines (with speed — (3 X while the remaining 14 propagate with non trivial speeds some of which depend on the gauge 
function /bm(oO- Again, it is important to emphasize that the strong hyperbolicity arises only with the value 0^0, 
otherwise it is not possible to construct all the eigenfields. Moreover, with the optimal value = 1, together with a 
positive definite NMC function /(</>), the eigenfields are all smooth in their variables. We remark = 1 is also the 
required optimal value in the second order approach (Cf. Ref. [16| ) for the reduced system of field equations to be in 
quasilincar diagonal second order hyperbolic form m = 1 was assumed in Ref. [la ] , and the shift was also included 
as a live shift). 

Basically both the second and the first order analysis sketched here show that STT posses a well posed Cauchy 
problem in the Jordan frame. 



IV. STATIC AND SPHERICALLY SYMMETRIC INITIAL DATA 



In this section we shall be concerned with initial data that represent a static and spherically symmetric (SSS) 
spacetime. In fact, we will be interested mainly in initial data corresponding to static but unstable configurations 
of compact objects (namely, neutron or boson stars). Here by unstable we mean a configuration that is near to 
the transition to spontaneous scalarization [2?], [H, [23, H(j or near to the point (i.e., maximal mass) of gravitational 
collapse towards a black hole. In this way, a small perturbation to the unstable configuration will lead to a non trivial 
dynamics that can be accompanied by the emission of scalar gravitational waves [29( . 

First, staticity implies that the extrinsic curvature and the shift vector will be globally null. Moreover, the field 
variables will be time- independent. In this case the l.h.s of the momentum constraint (|16[) becomes identically null 
which implies J] 113 -** = 0. The remaining field equations read as follows: 
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This slicing condition is a bit modification with respect to the one considered in [Ttl in that the term contains the factor /bm m 
the denominator. However, both conditions agree when adopting the modified harmonic condition (termed pseudo harmonic in [id ) 

e = i = / BM . 
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Moreover, using Eq. (f2"6"]) one can replace the terms with 13; in Eqs. (|2"3")l and (|23|) to obtain respectively 
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(28) 



where E = n a n b T a i, is the effective energy-density which incorporates all the energy-density contributions (including 
the scalar field and the NMC). For instance, in the minimal coupling case / = l/(87rGo) = const., therefore E — 

E ma tt + 2 — ^ ^(^) wri ich corresponds to the energy-density of an ordinary scalar field plus the energy density of 
matter. 

Equations ([28]) and ([26]) are elliptic equations for the lapse and the scalar field respectively (notice D\Q l = 3 A</>). 



Finally, the conservation equation 



(29) 



together with the staticity assumption will provide the equilibrium equation for the matter which we have not specified 
yet. 



A. Spherical symmetry 

One of the simplest applications of the field equations of STT is for SSS spacctimes, which can be described by the 
following metric in area-r coordinates 

ds 2 = -a 2 {r)dt 2 + A 2 (r)dr 2 + r 2 d6 2 + r 2 sin 2 (6>) dip 2 . (30) 
Now, adopting the parametrization A 2 (r) = f 1 — 2G o™( r ) ^ ; onc can easily verify that Eq. (|2"7| reduces to 

d r m = 4irr 2 E , (31) 

where E is to be read-off from Eq. ([27} with ^Dia = (Q r d r a)/(aA 2 ) and Q 2 = (Q r ) 2 /A 2 . The integration of Eq. 
([3"T]) provides the "mass function" m(r) which at infinity (and for asymptotically flat spacetimes) corresponds to the 
ADM-mass of the configuration. 
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For the lapse we can use Eq. (|28p. However, it turns that the 9 — 9 component of Eq. provides a first order 
equation for the lapse. In order to arrive to that equation we use 3 R e s = Gom/r 3 + (God r m)/r 2 = Gom/r 3 + 4irGoE. 
Then one can verify that Eqs. (pM|) and (f2"T)) together with the explicit expression for E lead to 
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Finally Eqs. (j26f and (fTTj) read respectively 
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(33) 



d r (j) = Qr 



(34) 



Equations (13~T]) — (f3"4")l together with Eq. (|2T?|) are the set of fundamental equations to be solved for a SSS spacetime 
in the area-r coordinates. Notice that when one replaces the term d r a/a appearing in E and in Eq. (|33f by the 
r.h.s of Eq. (13"2")) . the fundamental equations have all the form d r U = J$?(r, U) where U stands for the field variables 
(</>, Q r , hia, m), and Ji? are non linear functions of their arguments but independent of the derivatives of the fields U. 
One can therefore employ a Runge-Kutta algorithm to solve the equations numerically. As far as we are aware this 
system of equations has not been considered before, except for very particular cases. In fact for the specific class of 
STT where / = (1 + 16nGn£6 2 ) / (8nG ) and T£ fc att = {p + p)u a u b + g ab p is the EMT of a perfect fluid, Eqs. (|5T])-(|531) 
reduce to those of Ref. [28| where the phenomenon of spontaneous scalarization (2?| was further analyzed. For a 
perfect fluid Eq. (f2"5]) yields the equation of hydrostatic equilibrium 



d r P = ~(p + p) 



(35) 



Moreover, 5^ tt r 



p, -E ma tt = P, T ma ^t = 3p — p. This differential equation can also be written in terms of the 
baryon density n(r) since for neutron stars the equation of state is usually parametrized as p — p{n), p = p(n). We 
have 



d r n 



(p + p)n d r a 
Tp a 



(36) 



where T := din [p] /din [n] is the adiabatic index. Then n(r) adds to U as a further field. 

As one can easily observe, in the absence of a scalar field [i.e.,/ = l/(8irGo) — const., = = V((f>) = Q r ], Eqs. 
(l3"Tj) and (|3"2"|) with (|33|) or ([3^)1 are the equations of hydrostatic equilibrium for a relativistic star in pure GR. 



The system of equations presented in this section can be used to construct unstable SSS initial data (compact 
objects) that can serve to trigger the dynamical transition to spontaneous scalarization or to trigger the gravitational 
collapse of an scalarized object towards a black hole formation (cf. Ref. One expects that if the static 

configuration is very close to the point of instability, a small perturbation to the unstable static solution will drive the 
configuration to the final stage which can be a scalarized object (if the initial configuration is an unscalarized neutron 
star near to the point to the phase transition) or a black hole (if the initial data represent a scalarized neutron star 
with maximal mass). 

Another possibility is to use a boson field as a matter model [3(| ■ In that case, the equation of hydrostatic 
equilibrium (j35|) is replaced by the Klein-Gordon equation associated with the boson field. 



V. ISOTROPIC AND HOMOGENEOUS SPACETIMES 



Another simple application of the 3+1 system of equations of Sec. II is for a Friedmann- Robertson- Walker (FRW) 
cosmology where the metric is given by 



ds 2 



-dt 2 + a 2 (t) 



dr 2 
1 — kr 2 



+ r 2 d9 2 + r 2 sin 2 (9) dip" 



(37) 
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which describes an isotropic and homogeneous spacetime. In this case the extrinsic curvature is given by K l j — — 2 <$ 8 . ; 
and so K = —3- = —3H(t), where Hit) is the Hubble expansion rate, and 3 i? = 6k/a 2 (where k = ±1,0). Here a 
dot indicates derivative with respect to cosmic time t. The Hamiltonian constraint (| 1 5 j) reduces to the Friedmann 
equation 



H -\~ ^ = — — T 



n 2 

-Bmatt + V{<S>) + — - 3/'fflI 
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On the other hand Eq. (|18[) provides the dynamic equation that governs the acceleration of the Universe 
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Finally, Eqs. (fT2|) and (| 19[) take respectively the form 

4> = n , (40) 



Eqs. P8j) — (|4T|) together with H — a/a are the fundamental equations that determine the dynamics of a homoge- 
neous and isotropic Universe. Like in previous section, we emphasize that in the absence of a scalar field Eqs. (|38|) 
and (I39|) reduce to the very well known equations in GR. 

Assuming a perfect fluid for the matter we have again ^matt = 3p, -E ma tt = P, Tmatt = 3p — p. The conservation 
Eq. (JUD then yields 

p + 3(p + p)H = , (42) 

which for a non interacting two fluid component (photons and pressureless matter) provides the matter energy-density 
as a function of the scale factor: p — ci/a 4 + C2/a 3 where c\ and C2 are constants that are fixed, for instance, with 
the values of the photon and matter densities today. 

As in the previous section, the cosmological equations (13"9")) — (T4"T]) have all the form U = J%?(t, U) 7 where now U 
stands for the field variables (0, II, a, H), and Jif are non linear functions of their arguments but independent of the 
time derivatives of the fields U. The form of the equations again permits the use of a Runge-Kutta method to solve 
the system numerically. In this respect, when the system is evolved backwards in time by taking as initial data the 
conditions of the Universe today, the Hamiltonian constraint (|38|) can be used to check the numerical accuracy at 
every time step. In particular, the initial data cannot be arbitrary but has to be consistent with Eq. (|38p as well. 

Indeed for the particular class of STT given by / = (1 + 167rGo£</> 2 )/(87rGo) a similar system of equations was 
used in (3lj to analyze in detail several cosmological models which can be viewed as models of the Universe with an 
effective gravitational "constant" that evolves in cosmic time. 

For the case with an arbitrary /(</>), equations similar to the system ([38 ]) — (|41 |) have been analyzed in the literature 
[121 [Tij although not obtained from a 3+1 formalism and therefore they are presented in a more complicated fashion. 
Using the Brans-Dicke parametrization, an analogous system of equations can be seen in Ref. (l7j . 



Actually the cosmic time does not appear explicitly in the equations. However, it is often convenient to use r := ln[a(t)/ao] (where 
ao is the scale factor today) as "time" variable instead of t. Then the equations take the form U = Jff{r, U), where now U stands for 
(0, n, H ) and the dot indicates derivative with respect to r. In this case, the explicit dependence of J£°(t, U) with respect to r arises 
because the matter energy-density, p, depends explicitly on the scale factor a. 
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VI. OUTLOOK AND CONCLUSIONS 



We have presented the scalar tensor theories of gravity in the Jordan frame and wrote their corresponding 3+1 
equations (constraint equations and evolution equations) in order to analyze the Cauchy initial value problem. Then 
we briefly discussed the procedure that led to a set of (strongly) hyperbolic equations for which a well posed Cauchy 
problem can be established. This analysis shows that the STT in the Jordan frame have a well posed Cauchy problem, 
contrary to what it is usually claimed in the literature. Two simple applications of the equations were presented. 
First the case of static and spherically symmetric spacetimes where a set of "master" equations were obtained for 
arbitrary NMC functions f (</>)■ These equations are specially suited for a numerical analysis. For a specific class of 
STT such equations have in fact showed their usefulness [28|. This situation is very interesting for the study of the 
phenomenon of spontaneous scalarization that arises in compact objects. Moreover, configurations of this kind at 
the threshold of instability (towards spontaneous scalarization or towards the collapse into a black hole) can serve 
as a physical initial data for studying the dynamics of such objects in STT. For instance, recent studies show that 
spontaneous scalarization can have an impact on the signature of the wave forms emitted by a perturbed neutron star 
|32| . Such wave forms when compared with future observations can be useful to validate some classes of STT or to 
bound the non-minimal coupling constants. Finally, we discussed the case of isotropic and homogeneous spacetimes 
(FRW cosmology) in STT. We again obtained a set of fundamental equations which govern the dynamics of the 
Universe in the context of STT. These equations have been quite useful for the numerical analysis of cosmologies with 
specific choices of /(</>) [3l| . 

In the near future we plan to use the formalism presented here to study the dynamics and transition to spontaneous 
scalarization in different kinds of compact objects. Moreover, in the cosmological setting we expect to tackle the 
acceleration of the Universe and the dark energy problem without the introduction of a cosmological constant but 
in the framework of STT. This issue has been extensively analyzed in the past (see Refs. [lj], HH for a review), 
nevertheless we consider that the formalism presented here allows to appreciate the NMC contributions in a more 
transparent fashion. 
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APPENDIX A: SIMPLE EXAMPLE OF AN HYPERBOLIC SYSTEM 

Let us consider the system 

d t Tl + ad x n + bd x Q = S II {Q,TL) , (Al) 
d t Q + cd x Tl + dd x Q = S Q {Q,Tl) , (A2) 

where the coefficients a — d are in general functions of (£, x) and Su,q(Q, II) are source functions. 

Using the judicious guessing approach we shall construct a new system of the form (|22[) which is manifestly hyper- 
bolic and then we shall confront this method with the more systematic method which uses linear algebra. 

First we ask for a linear combination w := 11 + oQ to be the eigenfunction. This means that we write an evolution 
equation for w, where a(a, b, c, d) is a function of the coefficients a — d which is to be found such that the r.h.s of the 
evolving part take the form X a d x (II + o~Q), where X a will be the eigenvalue(s). In this way we have 



d t (n + o-Q) 



-d x (a + ac) II + (b + ad) Q 
b_ + ad Q 



{a + ac) d x n 



(A3) 



Note that with the symbol ~ we have discarded from the analysis all the terms which do not contain derivatives 
of II and Q (including the sources). In this sense we can consider as though the coefficients a — d were constants. 
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So, by comparing the coefficients of Q in the evolving <9 t and the spatial d x parts of Eq. (|A3|) . and then imposing 

er = (b+ ad) /{a + crc), one finds a± = - a± v / (°^ d ) + 4& ^ jy[ oreover ^ ^he eigenfunctions will be smooth provided c(t, x) 
does not vanish in some spacetime point. Therefore we conclude that the two eigenfunctions are w± = Tl + a±Q with 
the corresponding eigenvalues X± = a + a±c = a + rf± V(° d ) + 4b ^ rp^g e ig enva l U es are rea j arLC j the eigenfunctions are 

non degenerate if (a — d) 2 + 46c > and c 7^ 0. If this condition holds and the eigenfields are smooth then the system 
is strongly hyperbolic 8 . 

Now, let us find the eigenfunctions of the same system using the standard method of linear algebra. We rewrite 
the system in a more convenient matrix form as: 



The eigenvalues A of the characteristic matrix turn to be A± = a+rf± %/(° d ) + 4bc ^ w j 1 j c ] :l are identical to the ones found 
before. On the other hand, the eigenvectors are: 



u± = ^% . (A5) 



-1 

Note that the first component is identical to a^. Thus, the diagonalizing matrix reads 

<7_ (T + 
v -1 -1 

Then we can write the eigenfunctions as: 



h J \Q J a+ - a- \ n + a -Q 



(A6) 



(A7) 



which are the same found before using the judicious guessing approach modulo factors - T i- = T C /"\J (a — d) 2 + Abe. 
The fact that we do not obtain exactly the same eigenfunctions is because the eigenvectors are unique up to a 
rescaling. However, both couples of eigenfunctions w± and / 12 are physically equivalent. Finally, notice that if a = d 
then A± = a ± Vbc, and w± = II ± y^b/cQ thus 'a' will play the same roll as — (3 l in the field equations of Sec. 2. 
For instance, if the coefficients b, c are null, the equations (|A1|) and (|A2[) decouple and then n and Q are themselves 
eigenfunctions, indicating that they propagate along the "time lines" (with speed a). On the other hand, when b and 
c are non null, these coefficients are related to the propagation of the eigenfields along the "light cones" . For instance, 
b and c can be identified with a^f xx and a, respectively, in many of the field equations (assuming propagation only in 
the 'a;' direction). Thus \± are identified with the speeds of propagation — X ± a^/j xx (along the line cones) of many 
of the eigenfields used in the hyperbolicity analysis of STT [cf. Eq. (|22|) ] [25| or even GR Q (e.g. w± = u\ ± \J^f xx ui 
where u\ and represent two field variables of STT [cf. Eq. ([2H)l ]). 
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